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We consider particle-like and black holes solutions of the Einstein- Yang-Mills system with positive 
\ cosmological constant in d > 4 spacetime dimensions. These configurations are spherically symmetric 

. and present a cosmological horizon for a finite value of the radial coordinate, approaching asymptotically 

' the de Sitter background. In the usual Yang-Mills case we find that the mass of these solutions, evaluated 

outside the cosmological horizon at future/past infinity generically diverges for d > 4. Solutions with 
finite mass are found by adding to the action higher order gauge field terms belonging to the Yang-Mills 
hierarchy. A discussion of the main properties of these solutions and their differences from those to the 
usual Yang-Mills model, both in four and higher dimensions is presented. 

o 

O '. 1 Introduction 

cr 

, Recently, there has been a significant increase in interest in the properties of gravity in more than d = 4 

dimensions. This interest was enhanced with the development of string theories, along with the idea of 
large extra dimensions recently resurrected by TeV gravity models. Several solutions of higher dimensional 
classical general relativity have been known for some time, in particular extensions to any d > 4 of the 
Schwarzschild and Reissner-Nordstrom black holes by Tangherlini [T], and of the Kerr black hole by Myers 
and Perry [5] . Investigations over the past fifteen years have produced an impressive catalogue of solutions 
for various effective theories of Einstein gravity coupled to many different kinds of matter fields. These 
results indicate that the physics in higher-dimensional general relativity is far richer and complex than in 
the standard four-dimensional theory. 

Solutions to the Einstein- Yang-Mills (EYM) equations in higher dimensions have recently been studied. 
As found in [3J, for asymptotically flat solutions to the usual Yang-Mills (YM) gravitating system in five 
spacetime dimensions the particle spectrum obtained by uplifting the d = 4 flat space YM instantons become 
completely destroyed by gravity, as a result of their scaling behaviour. These results [3J can be systematically 
extended to the d > 5 case and one finds that no finite mass spherically symmetric solutions exist in EYM 
theory, unless one modifies the non Abelian action density by adding higher order curvature terms in the 
YM hierarchy. (The YM hierarchy features higher order curvature 2p forms F(2p) = F A F A F... A F, p 
times, labeled by p, the p = 1 term giving the usual YM system.) Without these higher order YM terms, 
only vortex- type finite energy solutions [3J exist, describing effective systems in three spacelike dimensions 
and with a number of codimensions. 

Asymptotically flat, regular, static and spherically symmetric solutions of EYM equations with higher 
order terms in the Yang-Mills hierarchy were presented in j¥j for spacetime dimensions d = 6, 7, 8, and for 
d = 5, both globally regular and black hole solutions were found in [5]. The properties of these solutions 
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are rather different from the familiar Bartnik-McKinnon solutions [6\ to EYM equations in d — 4, and 
are somewhat more akin to the gravitating monopole solutions to EYM-Higgs system [7]. This is because 
like in the latter case [7J, where the vacuum expectation value of the Higgs field features as an additional 
dimensional constant, here also additional dimensional constants enter with each higher order YM curvature 
term [3] . The typical critical features discovered in [U [5] have been analysed and explained in [5] . 

When a cosmological constant is added to the theory, the asymptotic behaviour of the physically relevant 
solutions changes from Minkowski to Anti-de Sitter (AdS) for A < or de Sitter (dS) for A > 0. In the 
latter case (A > 0), to our knowledge no studies of EYM solutions in higher dimensions (d > 4) have been 
undertaken to date. To carry out this investigation is the aim for the present work. 

Higher dimensional EYM-F(2) solutions with a negative cosmological constant (A < 0) have been studied 
in [21 [TO]- The main properties of these configurations resemble the familiar d = 4 AdS ones [TTJ IT2"] . 
These describe a continuum of solutions with arbitrary asymptotic values of the function w(r) parametrising 
the spherically symmetric gauge field. However, the total mass-energy of the AdS^ nonabelian solutions 
diverges for d > 4. Higher dimensional asymptotically AdS EYM solutions with finite mass-energy were 
found by augmenting the action density of the system with higher order curvature terms, consisting of 
2p-form curvatures F(2p) 10J. The qualitative features of these solutions are very similar to those of 
the asymptotically flat case with A = 0; in particular, and most notably, the nonabelian fields approach 
asymptotically a pure gauge configuration with w(r — > oo) = —1, uniquely, unlike in the p = 1 case. 

Proceeding to EYM systems with positive cosmological constant A > 0, we note that at present only 
the p = 1 gravitating YM system is studied, and that, only in d — 4 (see [T3]-[TB], and also the systematic 
approach in [17]). The gravitating YM field equations present both black holes and solutions with a regular 
origin, in contrast to the case of a gravitating Abelian field. This property is shared with all cases, A > 0, 
A < and A — 0. Asymptotically dS (A > 0) configurations present, in particular, a cosmological horizon for 
a finite value of the radial coordinate. As in the A < case, the asymptotic value of the gauge field function 
w(r) for solutions with A > is not fixed, implying the existence of a nonvanishing magnetic charge. This 
contrasts with the asymptotically flat situation [6]. 

Our task in this paper is to examine the corresponding situation for d > 4 EYM solutions with positive 
cosmological constant. Our strategy is to first consider the usual YM model, namely the p = 1 member of 
the YM hierarchy, i.e. the square of the 2-form curvature F(2). We find that although the EYM equations 
in this case present solutions approaching asymptotically the dS background, the mass of solutions evaluated 
at future/past infinity generically diverges. Like in the asymptotically flat and AdS cases, finite mass-energy 
EYM solutions are found by augmenting the action density of the system with higher order curvature terms, 
consisting of 2p-form curvatures F(2p). 

2 Higher dimensional gravitating p = 1 YM system 
2.1 The model 

In this Section we shall examine the usual EYM system in a d— dimensional spacetime described by the 
following action 

/ = hulk + Isurf = [ d d x^(-^—(R-2A) + C m )--^—[ d^xVhK, (1) 

where R is the Ricci scalar associated with the spacetime metric g u[/ , A — (d — l)(d — 2)/(2£ 2 ) is the 
cosmological constant and G is the gravitational constant (following [H[S], we define also k = l/(8irG)). 
dM ± are spatial Euclidean boundaries at Euclidean surfaces at future/past timelike infinity X ± and f gM ± 
indicates the sum of the integral over the early and late time boundaries. The quantities t^,„, h^ v and K are 
the bulk spacetime metric, induced boundary metrics and the trace of extrinsic curvatures of the boundaries 
respectively. 

The matter term in the above relation, C m = —\t\ tr {F^F^} is the usual F(2) nonabelian action 
density, F^ v = d^A v — d u A^ — ilA^, A v \ being the gauge field strength tensor. 
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The field equations are obtained by varying the action (Q} with respect to the field variables g^vjA 
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R„ v - \g„uR + A 5m „ = 8nGT^, £)„ (V=5 F"") = 0, (2) 

where the energy momentum tensor is defined by 

= tr {F^Fvpcf* - ^.F^F^}. (3) 

For the case of a d-dimensional spacetime, we restrict to static fields that are spherically symmetric in the 
d — 1 spacelike dimensions, with a metric ansatz in terms of Schwarzschild coordinates 

dv^ 

ds = ivw + r2 ^~ 2 " ° 2 ^ N ^ dt2 > W 

with dfld-2 the d — 2 dimensional angular volume element and 

the function m(r) being related to the local mass-energy density up to some d— dependent factor. 

The choice of the gauge group compatible with the symmetries of the line element ((4J is discussed in [4j . 
This choice implies the use of the representation matrices SO±(d), where d = d and d — d — 1 for even and 
odd d respectively. In this unified notation (for odd and even d), the spherically symmetric Ansatz for the 
S'0±((i)-valued gauge fields then reads [4] 

!- w ( r A v (±)^' v(±) _ _1 ( 1 ± F ^+i 



The T's denote the d-dimensional gamma matrices and 1, j = 1, 2, d — 1 for both cases; x 3 — x J /r, with 

2 ? 

r = Xja; . 

Inserting this ansatz into the action ((T|) , the EYM field equations reduce to 

(r d - A aNw')' - (d - 3)r d -<W - 1) W = 0, m' = ^r d ~ 4 f A^' 2 + (d - 3) ^ ~ , - = H— , (7) 
x ' 2 \ 2r^ / cr k r 

where the prime denotes the derivative with respect to the radial coordinate r 

The above differential equations have two analytic solutions. For a pure gauge field w(r) = ±1, one finds 
m(r) = M, o~(r) = 1, M being a constant, which corresponds to Schwarzschild-dS spacetime. For w(r) = 
we find a non Abelian generalisation of the magnetic- Reissner-Nordstr0m-dS (RNdS) solution with <r(r) = 1 
and 

m(r) = M + ?± log(J) for d = 5, and m(r) = M + r d ~ 5 for d ^ 5, (8) 

Mo being an arbitrary constant. Some properties of these solutions are discussed in [18]. We can see e.g. 
that, for a suitable range of (Mo, A) it describes a cosmological black hole, the horizons being located at the 
zeros of N. Also, although these solutions are asymptotically dS, their total masses/energies defined outside 
the cosmological horizon at future/past infinity (r — > oo), diverge. 

2.2 Numerical solutions 

We want the generic line element ^) to describe a nonsingular, asymptotically de Sitter spacetime outside 
a cosmological horizon located at r = r c > 0. Here N(r c ) = is only a coordinate singularity where all 
curvature invariants arc finite. 
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Figure 1. The functions cr(r), w(r), N(r) and m(r) are plotted as functions of radius for a typical one-node 
d — 5 regular solutions in a F(2) EYM-dS theory. 
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Figure 2. The value of the cosmological horizon radius r c , the value M in of the mass function m(r) at r — r c , 
the Hawking temperature associated with the cosmological horizon as well as the value a(0) of the metric function a 
at the origin, are shown as functions of A for d = 5 particle-like solutions of F(2) theory. 

Outside the cosmological horizon r and t changes the character (i.e. r becomes a timelike coordinate for 
r > r c ). A nonsingular extension across this null surface can be found just as at the event horizon of a black 
hole. The regularity assumption implies that all curvature invariants at r — r c are finite. Also, all matter 
functions and their first derivatives extend smoothly through the cosmological horizon, e.g. in a similar way 
as the U(l) electric potential of a RNdS solution. 
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Figure 3. The functions a(r), w(r), N(r) and m(r)/r are plotted as functions of radius for a typical d = 6 
regular solutions in a F(2) EYM-dS theory. 
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Figure 4. Typical EYM-f(2) solutions in d — 8 asymptotically flat and AdS spacetimes. 

As in the A < case, it is natural to consider two types of configurations, corresponding in the usual 
terminology to cosmological particle-like and black hole solutions. The black hole configurations possess an 
event horizon located at some intermediate value of the radial coordinate < < r c , all curvature invari- 
ants being finite asr-»r/,. Both the event and the cosmological horizons have their own surface gravity k 
given by 



K h,c 



-9 tt g rr {d r g tt f 
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Figure 5. The profiles of two different d — 5 EYM-f(2) black hole solution are presented for the same value 
of the cosmological horizon radius. 

the associated Hawking temperatures being T^' c = |k/ 1jC |/(27t). 

The corresponding boundary conditions at the origin and cosmological/event horizon are found by taking 
T2 = in the general relations (fl8|) . IT9|) given in Section 3. 

To integrate the field equations, we used in this work the differential equation solver COLSYS which 
involves a Newton- Raphson method [T§]. The case A > leads to the occurrence of a cosmological horizon 
at r = r c with N(r c ) = 0. As in the U(l) case, the cosmological horizon radius r c is a function of A. In 
practice, we solved the equations first on the interval [0, r c ] (or [r^,r c ]), choosing r c by hand and imposing 
regularity conditions of the solutions at r = r c . This allows to determine the values of the functions 
(m,a,w,w') at r = r c as well as the numerical value of A corresponding to the choice of cosmological 
horizon. The integration on [r c , oo) can then be performed as a second step, leading to the knowledge of the 
solution on the full r-axis and regular at r — r c . 

For a F(2) theory, the constants K and t\ can always be absorbed by rescaling r — ► cr, A — > A/c 2 and 
m — ► totic^ -5 , with c = y t~i/k, the only remaining parameter being the cosmological constant. We have 
numerically integrated Eqs. {7]) for d = 5,6, 7, 8, and for several values of A. 

All numerical solutions we found have w 2 (r) < 1, which implies a nonzero node number of the gauge 
potential w. This can be proven by using the sum rule 

- rf " 4 AW\ r c = J\ ra ^ d _ 3 y- 6{1 _ w 2 ) + r d-A N ^y (9) 



w 



which follows directly from the YM equations (ro here is ro = for solutions with regular origin or ro = , 
for black holes). Suppose that w(r) never vanishes and w 2 < 1 for ro < r < r c . Then the l.h.s. of the above 
relation vanishes, while the integrand of the r.h.s. is positive definite. Therefore the gauge potential of the 
nontrivial YM configurations with w 2 < 1 must vanish at least once in the region inside the cosmological 
horizon. For solutions with dS asymptotics, the asymptotic value of the gauge potential w(oo) = wq appears 
as a result of the numerical integration. As a general feature, all configurations we have considered have 
Wo < 1- 

Considering first the case d = 5, we were able to construct a numerical solution for each value of the 
cosmological constant A < 0.5, the numerics becoming more difficult for larger values of A. As in the d = 4 
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case [13], we expect there to be a maximal value A c of the cosmological constant, such that for A > A c the 
gravitational interaction becomes too strong for the solutions with dS asymptotics to exist. The profile of 
a typical one-node solution corresponding to A = 0.1 and r c = 3 is given on Figure 1 (multi-node solutions 
have been found as well). One can see that the gauge field function approaches asymptotically a finite value, 
Wq ~ —0.98 as a consequence of which the mass function diverges logarithmically. 

The corresponding numerical value of r c , u(0), the Hawking temperature Tjj of the cosmological horizon 
and the value of the mass function on the cosmological horizon M in — m(r c ) are reported on Figure 2 as 
functions of the cosmological constant for one- node configurations. One can see that, similar to the case of 
a RNdS solution, the Hawking temperature increases with A. 

Extending these solutions outside the cosmological horizon reveals, however, that the value wo differs 
slightly from one. As a general feature, the Einstein equations imply that the mass function m(r) develops 
a logarithmic term which makes the mass divergent as r — ► oo, 

m(r) = M + ^(wl-lf\og{j). (10) 

A similar result holds for d = 5 asymptotically AdS solutions [10] . the A = picture being more complex 

The results we found by solving numerically the field equations for d — 6, 7, 8 confirm that this is a 
generic behaviour of the higher dimensional EYM solutions. For d > 5 the mass function diverges for large 
r according to 

-M = M„ + ^f|K 2 -i)V- 5 , (ii) 

other properties of these solutions being very similar to the the d — 5 case. Despite this divergence, these 
solutions are still asymptotically dS (as r — ► oo one finds a ~ 1 + 0(l/r 6 ) ). A typical d — 6 configurations 
with a regular origin is presented in Figure 3, for A = 0.0001 and r c — 77.45. One can see that the mass 
function diverges linearly, while the gauge potential w(r) presents three nodes. 

For completeness, we give in Figure 4 the profiles of typical EYM-F(2) d = 8 configurations in asymptot- 
ically flat and AdS spacetimes (the solutions for other dimensions d > 5 have the same qualitative features). 
The solutions we have found for A < exist for a compact interval < b < b max , b being the parameter in the 
expansion of the gauge function at the origin w(r) = 1 — br 2 + 0(r 4 ), their masses diverging again according 
to l|lip. since Wp ^ 1. For the AdS case, the asymptotic value of the gauge potential may take arbitrary 
values [10], being fixed by the parameter b (with b = 0.5 for the solutions in Figure 4). All asymptotically 
flat solutions we have studied have wq = and present at least one node of the gauge function w(r). Also, 
a critical solution is approached as b — > b max , with er(0) — > in this limit. 

Apart from particle-like solutions, we have found black hole solutions as well. Restricting to the d = 5 
case, our numerical analysis suggests that any asymptotically dS particle-like solution presents black hole 
counterparts. Imposing the condition of a regular horizon at r — rh, we obtained a family of black hole 
solutions for any r c > rh > 0. When the value rh increases, the values |io(r/,)| and |w(r c )| slowly decrease. 
The asymptotics of the cosmological black hole solutions are similar to the particle-like case. In particular 
one finds Wq < 1, which implies a divergent value of the mass- function as r — > oo according to (jlOp . The 
profiles of the solutions corresponding two different values of are reported on Figure 5 (these configurations 
have the cosmological horizon at r c = 50). 

We see that the mass function of both regular and black hole solutions diverges asymptotically, yielding 
infinite total mass. The situation is exactly the same for the A < too. Not having a finite value for the 
mass in the A > case however may not be regarded as quite as serious a physical disadvantage, if one takes 
the view that the mass is nonetheless finite inside the cosmological horizon. 

For A < 0, the non-existence of d > 4 spherically symmetric EYM-F(2) configurations with finite mass 
could be proven in a rigurous way. However, the arguments in [9] [10] fail to apply for a positive cosmological 
constant, a different approach being necessary in this case. 
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3 Higher dimensional gravitating nonabelian p YM hierarchies 



3.1 The equations and boundary conditions 

A simple way to find nontrivial solutions with a finite mass is to modify the matter Lagrangean by adding 
higher order terms in the YM hierarchy, constructed exclusively from YM curvature 2p-forms. Such terms 
are also predicted by the low energy string theory (see e.g. [20]-|22j). 
The definition we use for superposed YM hierarchy is 



E 



-g Tr {F(2p) 2 } , 



(12) 



U 2(2p)l P 

where F(2p) is the 2p-form p-fold totally antisymmetrised product of the SO(d) YM curvature 2-form F{2) 

Even though the 2p-form ([15]) is dual to a total divergence, namely the divergence of the corresponding 
Chcrn-Simons form, the density (|12j) is never a total divergence since it is the square of one. But the 2p-form 
(I13p vanishes by (anti)symmetry for d < 2p so that the upper limit in the summation in (| 12[) is P = i for 
even d and P = for odd d. 



We define the p-stress tensor pertaining to each term in (fT2|) as 
T$ = Tr{F(2p) MAlA2 ... A 



1 F{2p) u x ^- x ^ - j-g^ F{2p) XlX2 .. Mp F{2p) x ^ - ^}. (14) 



We shall restrict in this work to solutions in dimensions less than nine, in which case it is sufficient to 
consider the first two terms in the YM hierarchy, i.e. a F(2) + F(4) model (see [HI [TO] for the equations of 
the general (P, d) model) . As in the previous section, we restrict our attention to static spherically solutions 
given by the Ansatze (J4| and ([6]), with exactly the same choices for the gauge group as in Section 2.1. 
The field equations of the F(2) + F(4) model are 



n ({r d - A CTNw')' - (d - 3)r d - 6 f j( W 2 - l)u>) + 
+H{d - 3)(d - 4)(w 2 - 1) ((r d - & crN(w 2 - l)w')' - (d - 5)r d - 10 a(w 2 - lfw^j = 0, 



(15) 



m! = -r d ^ | 



Nw' 2 + -(d-3) 



w 



+Ig(d-3)(d-4) 
6 r z 



Nw' 2 + -(d-5) 



w 



r 1 + l!|(d-3)(d-4) 



1 



w 



12 



The corresponding expansion of the gauge potential and metric functions as r — > is 

u>(r) = l-br 2 + 0(r 4 ), m{r) = (n + y (d - 3)(d - 4)6 2 )&V- 1 + 0(r d+1 ) 
2fo 2 (7o/ , 2t 2 , 



(16) 
(17) 

(18) 



cr(r) 



CO 



+ ^ (n + ^(d - 3)(d - 4)6 2 )r 2 + 0{r% 



and contains one essential parameter b (the value of <tq can be fixed by rescaling the time coordinate). 

Assuming the existence of a regular, nonextremal event horizon at r = ro (with ro = r/, or ro = r c ), the 
approximate expression of the solution near the event horizon is 
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m(r) = (1 - r d - 3 + m'(r )(r - r ) + 0(r - r ) 2 , (19) 

cr(r) = a Q + <r' (r - r ) + 0(r - r ) 2 , w(r) = u> + u/(r )(r - r ) + 0(r - r ) 2 , 

where 

■' 6 ,K 2 -1) 2 \ d-3 (d-l)ro 2m'(r ) 



m'(ro) = ^-(d- 3)(w 2 - l) 2 ( n + r 2 (d - 4)(d - 5) ^|^ ), ^ = 



^ 2 k^" 3 



o 7 r o «* KrJ 

°o = T i + r 2(d - 3)(d - 4) —. , w = — J (d- 3)( 2 " , ), 

V o o rg A Tl+T2 ( d _ 3 )( d _4)M^)! ; 

with two free parameters, wq, a$. Also, since the field equations are invariant under w — * —w, one can take 
w(0) = 1 and w(rh) > without any loss of generality. 

For r — > oo we find for both regular and black hole solutions 

.- M — -Li I w l ■ m(r) — M n(d-3)w 2 1 ^?(rf-3) 2 n l . 

w(rj -±1 + -^ + ..., m(rj-A7 ^ + . . . , cr(r) - 1 2K (d-2) + • • • • 

These boundary conditions are also shared by the asymptotically flat solutions (with a different decay of 
the mass function m(r), however), w = ±1 being again the only allowed values of the gauge function as 
r — > oo. As a general feature, all solutions discussed in the rest of this section present only one node in the 
gauge function w(r), i.e. w(oo) = —I. As in the A < cases, we could not find multi-node solutions in the 
F(2) + F(4) model. 

3.2 Numerical solutions 

In the presence of higher oder terms in the YM action, dimensionless quantities are obtained by rescaling 

r - (Ira/nl) 1 /^, A - (In/r^A, m(r) - K(\n/T 2 \)«- 3 V*m(r). (20) 

This reveals the existence of one fundamental parameter which gives the strength of the gravitational in- 
teraction a 2 = |ti| 3 / 2 /(k|t2| 1 / 2 ). The solutions can then be constructed in terms of a 2 and A. Most of 
the numerical work was carried out for positive coupling constants in the YM hierarchy (|12p . having set 
n = T2 = 1 without loss of generalitj0. After presenting our main results however, we briefly consider at 
the end of this section, the case where r 2 is taken to be negative. 

Starting again with particle-like solutions in d = 5, we have solved equations (|16p for several values of A, 
varying a. The pattern of these solutions is illustrated by Figure 6a, where the quantity <r(0) is reported as 
a function of a 2 for three different values of A. 

Here it is useful to recall the situation corresponding to A = (see [5]). In this case, several branches 
of solutions exist, depending on the parameter a 2 , as illustrated on Figure 6a. The first (or main) branch 
exists for a 2 £ [0,0.2824]. Then another branch of solutions (with larger masses than the corresponding one 
on the main branch) exists for a 2 E [0.1749,0.2824]. Several other branches of solutions further exist on 
smaller intervals centered on the critical value a 2 r ~ 0.1749. It is clear from Figure 6a that this oscillatory 
behaviour converging on a 2 r is a common feature of A > solutions. (This is also the case with A < 
solutions [10) . which is not displayed on Figure 6a since its branch structure is nearly identical to that of 
A = 0.) This phenomenon, which occurs in appropriately similar models in all 4p + 1 dimensions, was 
exhaustively analysed for the A = case in [8], where this a 2 T was called a conical critical point. In the 
present paper, we limit ourselves to a qualitative discussion only. 

For positive values of the cosmological constant, our numerical analysis in the case < A < 0.001 reveals 
that an extra branch of solutions exists, as shown on Figure 6a for A = 0.0004. 



1 The coupling constant ti equals the inverse of the square of the gauge coupling constant and is strictly positive. Also, 
string theory predicts a positive value for the coefficient of the -F(4) term in the YM hierarchy (see [20]- [22]). 
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Figure 6. Several global quantities are shown as functions of a 2 and several values of A for d = 5, 8 solutions 
of F(2) + F(4) theory. 

Our numerical analysis suggests that this branch, which has no counterpart in the A = case, exists 
for a > Oq, with w 0.21 for A = 0.0004. This branch, appears to survive for arbitrarly large values of a 
(although the numerical accuracy deteriorates for a larger than one). Decreasing the cosmological constant, 
the critical value ag decreases and the pattern of the A = case is approached. For larger values of A the 
pattern simplifies and a single branch persists, as illustrated for A = 0.0025 on Figure 6a, for a 2 > 0.1749. 
Also, only a small variation of the asymptotic value M of the mass function was noticed when varying the 
parameter a. 

Three remarks can be made to summarise this description : (i) For d = 5, A > solutions do not exist 
for arbitrarily small values of a; (ii) Non trivial cosmological solutions seem to exist for large values of a 2 , 
and (iii) The multiple branche phenomenon converging on the conical fixed point a 2 ~ 0.1749, observed for 
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Figure 7. The event horizon Hawking temperature, the value of the gauge potential on the event horizon Wh 
as well as the value of the metric function a at the event horizon, are shown as functions of the event horizon radius 
for d = 5 black hole solutions of F(2) + F(4) theory with a 2 = 0.25, A = 0.0004. 
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Figure 8. The functions cr(r), w(r), N(r) and m(r) are plotted as functions of radius for a typical d — 5 black hole 
solution in a F(2) + F(4) EYM-dS theory. 

A = 0, seems to persist for A > 0. 

The F(4) term allows also for the existence of d = 6, 7, 8 configurations with finite mass. The solutions 
available in the absence of a cosmological constant (A = 0) feature two branches which exist for a 2 £ [0, a£j, 
where the maximal value a m depends on d (see [4]). Integrating the equations for A > reveals that the 
cosmological solutions also feature two branches leading to a pattern very similar to the A = case. The 
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maximal values a m decrease with increasing cosmological constant. Also, with increasing a the masses of 
the gravitating solutions decrease. Along the second branch the value of a(0) decreases monotonically with 
a. The mass of a second branch solution is always larger than the corresponding mass (for the same value 
of a) on the first branch. Some relevant quantities for d = 8 solutions are plotted in Figure 6b. 

So far we discussed solutions which are regular at the origin. However the equations also admit black 
hole solutions with a regular event horizon occuring at with < r c . The study of the domain of black 
hole solutions in the space of the parameters (a, A, 77,) is a considerable task which is beyond the scope of 
this paper. We present however a few features of these solutions which reflect the general pattern, limiting 
our investigation to the d = 5 case. 

As we have seen in the previous section, several branches of regular solutions exist, according to the value 
of A. When imposing a regular horizon at r — r^, our numerical analysis indicates that the regular solutions 
are deformed into black hole solutions. We have analysed in detail the evolution of the black holes solutions 
in the case a 2 = 0.25, A = 0.0004 (corresponding to r c ~ 50.0). In this case there are three different regular 
solutions (see Figure 6) distinguished namely by the value of the metric function a at the origin (for instance 
<r(0) « 0.93, 0.79, 0.30, let us call them a, b and c respectively). The evolution of these solutions into black 
holes is summarized on Figure 7. We see clearly that the solutions b and c are deformed into black holes up 
to a rather small value of rh, namely for < 0.227. In the limit — > 0.227, the two branches merge into 
a single solution. A typical profile of a black hole solution corresponding to the solution c is presented on 
Figure 8. 

The scenario is completely different for the solution a. Indeed, it seems that it can be deformed into 
a black hole with large event horizon. When the value r/, increases, we observe that the corresponding 
function w(r) has a tendency to spead over the interval r £ [r^,r c ] and that the combination of w(r),w' (r) 
appearing in the right hand side of the equation for m'(r) becomes uniformly very small. As a consequence, 
the equation for m(r) can be simplified and leads to the following approximate solution 



for the region between the event and cosmological horizons. The metric functions N(r), <r(r) above turn 
out to be a very good approximations of the numerical solution a obtained for > 1 (although w(r) is 
non constant). Due to the presence of this solution, the numerical integration of the equation with the non 
trivial w(r) becomes increasingly difficult while increasing r^. 

The Hawking temperature associated with the cosmological horizon is almost constant on the branch (a) , 
and stronly decreases with for the other two branches. As a generic feature, the Hawking temperatures 
of the event and cosmological horizons are different (this holds also for the black hole solutions of the F(2)- 
theory). Therefore, the energy flows from the hotter horizon to the cooler one and the black hole will gain 
or lose mass. 

Before concluding this section, we allude to the unusual situation where the coupling constant T2 asso- 
ciated with the F(A) term is taken to be negative 0. We have managed to construct such solutions, which 
satisfy the same set of boundary conditions as the usual ones above with T2 > 0. Physically, this would 
imply a negative contribution to the mass-energy density of the solutions, but in practice this possibility 
does not seem to obtain and the solutions we constructed lead to finite masses. 

Since we do not have any existence proof for the type of EYM solutions found in this section, we are 
unable to exclude the possibility of solutions with T2 < 0. For this reason we have extended our numerical 
analysis to T2 < models with A < as well as the A > model at hand. Our numerical results indicate 
the absence of such configurations approaching asymptotically the AdS or flat background in all dimensions 
between five and eight. Surprisingly enough, we found finite mass solutions with T2 < in the dS case 
studied here. The results for d = 5 indicate the existence of a branch of solutions emerging at a — from 
a configuration in fixed dS background and extending in a up to a maximal value a ma x- This behaviour 

2 Ascribing a negative value to the square of the gauge coupling constant has been considered in the case of gravitating 
electromagnetism in [29) . where explicit solutions are found. We found however that in the corresponding situation when the 
Abelian field in 29 is replaced by a non-Abelian SU (2) gauge field, no numerical solutions with right asymptotics appear to 
exist. 





(21) 
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strongly constrasts with what we have found for r 2 > 0. However, the profile of a typical F(2) + F(4) 
solution does not depend on the sign of r 2 . As a — > a max the numerics deteriorates and the solver fails to 
converge (although all metric and matter functions stay finite there), a different approach being necessary. 
Along this branch, the value at the origin of the metric function a decreases and also the mass-parameter 
M (note that all solutions we found with r 2 < have M > 0, with a small variation when increasing a, 
however). The corresponding picture for r = 1, r 2 = — 1 and A = 0.0075 is plotted in the inlet of Figure 6a. 



4 A computation of the mass in the boundary counterterm 
method 

In evaluating expressions like the mass and action, one usually encounters divergencies coming from inte- 
gration over the infinite volume of spacetime. In the case of AdS gravity, a regularization procedure was 
proposed in [23] . that consists in adding to |T]) counterterms constructed from local curvature invariants of 
the boundary. These counterterms, which are esentially unique, can be easily generalized to the case of a 
positive A |24j . The following counterterms are sufficient to cancel divergences in a pure dS gravity theory 
for d < 7 



u- 1 



8irG 



d d ~h 



d-2 , m {d - 4) g _ fiejd- 6) _ d-i 2 



l 2(d-3) 2(d- 5)(d- 3) 2 V ^ 4(d-2) 



(22) 



Here IZab, TZ- are the Ricci tensor and the Ricci scalar for the boundary metric and Q(x) is the step function, 
which is equal to 1 for x > and zero otherwise; A, B, . . . indicate the intrinsic coordinates of the boundary. 
Jqm± indicates the sum of the integral over the early and late time boundaries. In what follows, to simplify 
the picture, we will consider the X + boundary only, dropping the indices (similar results hold for 

Using / = Ibulk + I surf + let, one can construct a boundary stress tensor, which is given by the variation 
of the total action at the boundary with respect to Hab (its explicit expression is given in Ref. |24|). 

The next step is to write the boundary metric in a ADM-like general form 

ds 2 = h AB dx A dx B = N 2 dt 2 + a ab [dijj a + N a dt) (dip b + N b dt) , (23) 

where N t and N a are the lapse function and the shift vector respectively and the ip a are the intrinsic 
coordinates on the closed surfaces E (a d— 2 dimensional sphere in our case). In this approach, the conserved 
quantity associated with a Killing vector on the I + boundary is given by 

0* = / d d - 2 i,^n A T AB e, (24) 



s 

where n A is an outward-pointing unit vector, normal to surfaces of constant r. Physically, this means that 
a collection of observers, on the hypersurface with the induced metric Habi would all measure the same 
value of 0£ provided this surface has an isometry generated by £ l . If d/dt is a Killing vector on E, then the 
conserved mass is defined to be the conserved quantity Tl associated with it. 

We have applied this approach to compute at the far future boundary (outside the cosmological horizon) 
the mass of the solutions of the F(2) + F(4) model. The crucial point here is that these solutions approaches 
asymptotically a Schwarzschild-dS background, the YM manifesting only in the next to leading order of the 
Tab expression. As a result, one finds 



(d-2)fi d _ 



SM =- - g"^'" 2 M + E (d), (25) 

where M is asymptotic value of the mass function m(r) and fld-2 = 2n( d ~ 1 ^ 2 /T((d — l)/2) is the area of 
a unit (d — 2)-dimensional sphere. The additional term Eo(d) appearing in (|25|) for d = 5, (7) is the mass 
of pure global dS 5 (7) spacetime and is usually interpreted as the energy dual to the Casimir energy of the 
CFT defined on a four (six) dimensional Euclidean Einstein universe [23] (i.e. E — 3£l 3 £ 2 / (64:irG) for d = 5 
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and Eq — Bil^l 4 / (128irG) in the seven dimensional case). One should also remark that all solutions of the 
F(2) + F(A) model we found have M > (both black holes and particle like solutions). Thus Wl — Eo(d) is 
negative, consistent with the expectation [25 that pure dS spacetime has the largest mass for a singularity- 
free spacetime. 

However, one can easily see that this approach fails to assign a finite mass to the solution of the EYM- 
F(2) model, despite their asymptotically dS behaviour. (Note that this prescription regularizes the mass 
and action of the embedded abelian solutions [53]). Asymptotically AdS solutions with a diverging ADM 
mass have been considered by many authors, mainly for a scalar field in the bulk (see e.g. 27 \). In this case 
it might be possible to relax the standard asymptotic conditions without loosing the original symmetries, 
but modifying the charges in order to take into account the presence of matter fields. For A < 0, the Ref. 
[TU] suggested that it is still possible to obtain a finite mass of EYM solutions in a F(2) theory by allowing 
I c t to depend not only on the boundary metric h A B, but also on the gauge field strength tensor. This means 
that the quasilocal stress-energy tensor also acquires a contribution coming from the matter fields. 

A similar approach holds also for dS solutions in d > 4 dimensions and we find that by adding to the 
expression (|22[) a supplementary matter counterterm of the form 

4 ( t m) = -T / d'xVhlogCj) tr{F AB F AB }, (26) 

4 JdM f- 

for d = 5 and 

= -nr^ I dd ' lxVK ^ f abF ab }, (27) 

l« _ °j JdM 

for d > 5, the mass divergence disappears. This yields a supplementary contribution to the boundary 
stress-tensor 

Tj$ = - T -^p- h AB tx{F CD F CD }, if d = 5, and = -J-J-h AB tx{F CD F CD }, for d > 5. 

The mass of the d > 5 solutions computed in this way is finite 

Tl=- {d ~ 8 2 ^ d ' 2 M Q + E (d), (28) 
where Mq is the constant appearing in the asymptotic expansion (jlip . 



5 Conclusions 

This work was primarily motivated by the question of how a positive cosmological constant will affect the 
properties of the gravitating nonabelian field solutions in a higher dimensional spacetime. To the best of 
our knowledge, this question has not yet been addressed in the literature. Apart from this motivation, the 
study of gravitating matter field configurations in asymptotically dS space may help a better understanding 
of the conjectured dS/CFT correspondence as well as act as a probe for the brane- world scenario. 

Our findings have completed our qualitative understanding of gravitating nonabelian solutions in higher 
(d > 5) dimensions, encompassing all possible values of the cosmological constant A < 0, A = and A > 0, 
the last being the object of the present investigation. Certain features of these solutions are shared, while 
others differ. At the most basic level, we have seen to varying degrees of rigour, that spherically symmetric 
solutions to the usual (p = 1), F(2) YM model in all these three cases have infinite mass in higher (d > 5) 
dimensions. By contrast, finite masses are obtained when the YM sector of the model is augmented by the 
appropriate higher order (p > 2), F(2p) members of the YM hierarchy in all three cases. (We have also 
considered the alternative option of using the counterterm method, avoiding the use of p > 2, F(2p) terms.) 
It can therefore be stated that, to construct finite mass solutions of gravitating nonabelian matter, the YM 
sector of the theory must be an appropriate superposition of members of the YM hierarchy, beyond the usual 
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YM term. This statement can be qualified in the present context, namely when A > 0, by adding that even 
when the total mass diverges the mass inside the cosmological horizon is finite. In addition to this salient 
feature of higher dimensional EYM solutions, our results reveal detailed qualitative properties occurring in 
all three cases. 

In the context of the p = 1, F(2) YM model with A > in d > 5 dimensional spacetime, the most 
remarkable property is that the asymptotics of the solutions lead to monopole-like configurations with 
nonvanishing magnetic flux. Indeed in those cases the solutions never have instanton like asymptotics. This 
circumstance makes it very easy to conclude that the total mass of these solutions is divergent. 

In the context of the F(2) + F(A) model studied in this paper, which is the simplest case of models with 
superposed F(2p) terms, our results have led to an overview of the main qualitative features common to 
EYM solutions in all three cases (A < 0, A = 0, A > 0), apart from their masses being finite. These confirm 
and expand on our knowledge of the branch patterns of these solutions in various dimensions, learnt from 
the results of 4,5 for A = 0, and [TU] for A < 0. The results of A = solutions were analytically analysed 
in [5] leading to a patterns that repeat modulo d = Ap + 1 dimensions. In particular there arise two types 
of patterns, those in dimensions d = Ap + 1 and the rest in Ap + 2 < d < Ap + A. In the restricted context 
here, these are the dimensions d = 5 on the one hand, and 6 < d < 8 on the other. We have learnt here that 
these patterns arise also for the A > case. In particular for the second case, solutions in d — 8 conform to 
the pattern displayed on Figure 6b in both the other cases. Much more interestingly the situation in d = 5, 
which is qualitatively the same for A < and A = models, strongly departs for A > 0, from the patterns 
of the former. These features are displayed on Figure 6a, and described in Section 3.2. 

What has not been studied quantitatively here, and in the A = [H [8] and A < [10] cases is the 
question of the stability of the solutions. We expect that in all even spacetime dimensions, as well as all 
odd dimensions d = Ap + 1, the solutions will be sphalerons like the four dimensional Bartnik-McKinnon 
solutions [5] . This is because in all these cases there is no topologically stable soliton in the gravity decoupling 
limit. In all other odd spacetime dimensions however, a stable soliton will survaive in the gravity decoupling 
limit, stabilised by the =^-th Chern-Pontryagin (CP) charge, rather like the monopole charge in the case of 
the gravitating monopole [3D], stabilised by the monopole charge (that descends from the 2-nd CP charge). 

As an additional remark, we allude to the case when the coupling strength of the F(A) term, T2, takes 
on a negative value. In the absence of existence proofs for solutions to models with higher order YM terms, 
one cannot a priori exclude this possibility, especially in view of the discovery of such explicit solutions in 
the usual Einstein-Maxwell theory [29] . We have tried to construct such solutions numerically and failed for 
models with A < 0, but surprisingly find them for the A > case at hand. 

While the present work is concerned with higher dimensional EYM solutions with positive cosmological 
constant, we have at every stage compared our results to the corresponding ones pertaining to the asymp- 
totically flat [H [5] and the asymptotically AdS [TU] counterparts. The general characteristics of both A < 
solutions, at least of the important finite mass ones, are quantitatively similar, while the corresponding 
characteristics of A > case differ from the former strikingly, in a systematic way. The numerically discov- 
ered features for A = solutions are explained analytically in [SJ, which can be adapted systematically to 
the A < case which is qualitatively similar. But an analytic study like [5] for the A > case, using the 
methods of [35J is outstanding and is desirable to complete the overall comparative study of all three types 
(A > 0, A = and A < 0) of higher dimensional EYM solutions. 
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